Nonanalytic behavior of the spin susceptibility in clean Fermi systems 
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The wavevector and temperature dependent static spin susceptibility, Xs(Qj T), of clean interacting 
Fermi systems is considered in dimensions 1 < d < 3. We show that at zero temperature Xs is a 
nonanalytic function of |Q|, with the leading nonanalyticity being |Q| d_1 for 1 < d < 3, and 
Q 2 ln|Q| for d — 3. For the homogeneous spin susceptibility we find a nonanalytic temperature 
dependence for 1 < d < 3. We give qualitative mode-mode coupling arguments to that effect, 

and corroborate these arguments by a perturbative calculation to second order in the electron- 
electron interaction amplitude. The implications of this, in particular for itinerant ferromagnetism, 
are discussed. We also point out the relation between our findings and established perturbative 
results for 1-d systems, as well as for the temperature dependence of X«(Q = 0) in d — 3. 
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I. INTRODUCTION 

It is well known that in fluids, that is in interacting 
many-body systems, there are long-range correlations be- 
tween the particles. For example, in classical fluids in 
thermal equilibrium there are dynamical long-range cor- 
relations that manifest themselves as long-time tails, or 
power-law decay of jeja.uilibrium time correlation func- 
tions at large times.Era In frequency space, the analo- 
gous effects are nonanalyticities at zero frequency. In an 
intuitive physical picture, these correlations can be un- 
derstood as memory effects: The particles "remember" 
previous collisions, and therefore so-called ring collision 
events, where after a collision the two involved parti- 
cles move away and later recollide, play a special role 
for the dynamics of the fluid. Technically, the long-time 
tails can be described in terms of mode-mode coupling 
theories. The salient point is that with any quantities 
whose correlations constitute soft, or gapless, modes (due 
to conservation laws, or for other reasons^, products of 
these quantities have the same property.u In the equa- 
tions of motion that govern the behavior of time corre- 
lation functions this leads to convolutions of soft prop- 
agators, which in turn results in nonanalytic frequency 
dependences. For phase space reasons, the strength of the 
effect increases with decreasing dimensionality: While in 
three-dimensional (3-d) classical fluids the long-time tails 
provide just a correction to the asymptotic hydrodynamic 
description of the system, in 2-d flwds they are strong 
enough to destroy hydrodynamics Iru 

A natural question to ask is whether such long-range 
correlations also occur in position space. Indeed, in clas- 
sical fluids in nonequilibrium steady states effects occur 



that may be considered as the spatial analogs of long- 
time-tails, but in thermal equilibrium this is not the 
case.oB This changes, however, if we consider quantum 
fluids. The quantum nature of a system has two major 
implications as far as statistical mechanics is concerned. 
First, temperature enters, apart from occupation num- 
bers, through Matsubara frequencies, which means that 
the system's behavior as a function of temperature will 
in general be the same as its behavior as a function of 
frequency, at least at asymptotically low temperatures. 
Second, and more importantly, in quantum statistical 
mechanics statics and dynamics are coupled and need 
to be considered together. This raises the question of 
whether in a quantum fluid there might be long-range 
spatial correlations even in equilibrium. 

From studies of systems with quenched disorder, there 
is evidence for the answer to this question being affir- 
mative. Let us consider interacting fermions in an envi- 
ronment of static scatterers. In dimensions d > 2, and 
for a sufficiently small scatterer density, the relevant soft 
modes in such a system are diffusive, so frequency f2, 
or temperature T, scales like the square of the wavevec- 
tor Q, ft ~ T ~ Q 2 . Via mode-mode coupling effects 
that are analogous to those present in classical fluids, 
dynamical long-range correlations lead to long-time tails 
in equilibrium time correlation functions. For instance, 
the electrical conductivity as a function of frequency be-i 
haves like £l (d ~ 2 ^ 2 at small Q, in d > 2 dimensions.tl 
The dynamical spin susceptibility, Xs(Q)^)i shows no 
analogous long-time tail at Q = for reasons related to 
spin conservation. However, from the above arguments 
about the coupling of statics and dynamics in quantum 
statistical mechanics and the scaling of frequency with 
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wavenumber, one would expect the static spin suscepti- 
bility, Xs(Q, H = 0) at T — 0, to show a related nonan- 
alyticity at Q = 0, namely \s ~ |Q| d ~ 2 - This is indeed 
the case, as,-can be seen most easily from perturbativc 
calculations Schematically, the coupling of two diffu- 
sive modes leads to contributions to Xs of the type, 
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which leads to the above behavior. One can then invoke 
renormalization group arguments to show that this is in- 
deed the leading small- Q behavior of \s- Similarly, at fi- 
nite temperature the homogeneous susceptibility behaves 
like Xs(Q = 0,Q = 0) ~ T (d - 2 ^ 2 . This has interesting 
consequences for itinerant magnetism in such systems, as 
has been recently discussedBtrB 

Somewhat surprisingly, the situation is much less clear 
in clean Fermi systems. Here the soft modes are den- 
sity and spin density fluctuations, as well as more gen- 
eral particle- hole excitations. All of these have a linear 
dispersion relation, i.e. SI ~ |Q|. The form of the disper- 
sion relation does not affect the basic physical arguments 
for nonanalytic frequency and wavenumber dependences 
given above. One might thus expect the spin susceptibil- 
ity to have mode-mode coupling cont ributions of a type 
analogous to those shown in Eq. (IT), but with ballistic 
instead of diffusive modes: 
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which leads to Xs(Q> ^ = 0) ~ const + |Q|" _i m generic 
dimensions at T = 0. In d = 3, one would expect a 
Q 2 In |Q| behavior, as convolution integrals tend to yield 
logarithms in special dimensions. Such a behavior of Xs 
would have profound consequences for the critical be- 
havior oLitinerant ferromagnets, as has been pointed out 
recently^ It is therefore of importance to unambigiously 
determine whether or not the above mode-mode coupling 
arguments do indeed carry over from disordered to clean 
systems. 

Before we start this task, let us discuss the available 
information concerning long-range correlations in clean 
Fermi systems. The specific heat is known to be a non- 
analytic function of temperature, viz. Cy/T ~ T 2 lnT 
in d = 3. This is a consequence of a nonanalytic correc- 
tion to the linear dispersion relation of the quasiparticles 
in Fefimi liquid theory, namely Ae(p) ~ (p — pf) In \p — 
pp\S3 Such a nonanalyticity signalizes the presence of a 
long-range effective interaction between the quasiparti- 
cles, and in general it will lead to nonanalytic behavior 
of both thermodynamic quantities and time correlation 
functions. The T 2 lnT term in the specific heat coef- 
ficient is an example of-such an effect. In d = 2 the 
behavior is Cy/T ~ TjTJ which is consistent with the 
behavior Cy/T ~ y<*— l j n g ener j c dimensions that one 
would expect from the above arguments. It was natural 
to look for similar effects in other quantities, in particular 



in the spin susceptibility. These investigations concen- 
trated on the temperature dependence of Xs , and several 
authors indeed reported to have found a T 2 In T term 
in the homogeneous static spin susceptibility. Howevec, 
other investigations did not find such a contribution.^ 2 ] 
The resulting-confusion has been discussed by Carneiro 
and Pethick.EJ These authors used Fermi liquid theory 
to show that, while T 2 lnT terms do indeed appear in 
intermediate stages of the calculation of Xs as well as of 
Cy, they cancel in the former. 

This somewhat surprising result casts some doubt on 
the general physical picture painted above, which sug- 
gests the qualitative equivalence of disordered and clean 
systems with respect to the presence of long-range corre- 
lations, and resulting nonanalyticities in both the statics 
and the dynamics of quantum systems. On the other 
hand, a failure of this general picture would be hard to 
understand from several points of view. For instance, in 
d = 1 the instability of the Fermi liquid with respect 
to the Luttinger liquid is well known to manifest itself 
in perturbation theory for Xs by means of logarithmic 
singularities .E5li3 This is precisely what one obtains from 
the mode-mode coupling integral, Eq. ( |l.2| ). By conti- 
nuity one therefore expects Xs(Q = 0,T) ~ T d_1 , and 
Xs{Q,T = 0) ~ |Q| d -\ at least in d = 1 + e. Unless 
the physics changes qualitatively between d = 1 + e and 
d = 3, this should still be true in higher dimensions. 
Also, the corrections to Landau theory we are discussing 
here can be cast in the language of the renormalization 
group. In this framework, the FerniiJiquid ground state 
is described as a stable fixed point Ji3 and the effects we 
are interested in manifest themselves as an irrelevant op- 
erator_that leads to corrections to scaling near this fixed 
point O In a system where Q, fl, and T all have a scale di- 
mension of unity, this operator should appear as |Q| d_1 , 
SI x , etc., dependences in various correlation functions. 
From a general scaling point of view it would be hard to 
understand if this were not the case, except for the possi- 
bility that the prefactors of some nonanalyticities might 
accidentally vanish in certain dimensions. 

It is the purpose of the present paper to clarify this 
confusing point. We will show that the above general 
physical picture does indeed hold true, and that it is 
not violated by the previously found absence of a T 2 In T 
term in Xs in d = 3, which is accidental. The remainder 
of this paper is organized as follows. In Sec. |n| we define 
our model. In Sec. [II we perform an explicit perturbative 



calculation to second order in the electron-electron inter- 
action. This confirms both our qualitative arguments, 
and the results of Ref. |l3|. We explain why there is no 
contradiction between these results, and we also make 
contact with established perturbative results in d = 1 . In 
Sec. IV A we discuss our result in the light of mode-mode 
coupling argum ents t hat are an elaboration of those given 
above. In Sec. IV B we make contact with renormaliza- 



tion group ideas, and argue that the functional forms of 
the nonanalyticities derived in Sec. [II by means of per- 
turbation theory are asymptotically exact. In Sec. IV C 
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we discuss the physical consequences of our results. 

II. MODEL, AND THEORETICAL FRAMEWORK 
A. The Model 



Let us consider a system of clean fermions governed by 
an actior£3 



S 



d 



J dxy^ y ip cr (x) ip.rir) + ,S', i + S ull . ! 2.1a.) 



Here we use a four- vector notation, x = (x, t), and 

f dx = J (ix J Q dr. x denotes position, r imaginary time, 
[3 = 1/T, and we choose units such that K = Ub = 1- cr is 
the spin label. So describes free fermions with chemical 
potential /i, 

S = f dxY^4><r{x) [A/2m + (i] ip a {x) , (2.1b) 

with A the Laplace operator, and m the fermion mass. 
Slut describes a two-particle, spin independent interac- 
tion, 



S lnt = ~ J dxxdxi ^2 v ( x i - x i) 

<J\,a-2 

X-0CT! (£l)VV 2 (X2)i><7 2 (^VVj (^l) 



(2.1c) 



The interaction potential v(x) will be specified in Sec. 
lTb| below. 



We now Fourier transform to wave vectors k and 
fcrmionic Matsubara frequencies tu n — 2nT(n + 1/2). 
Later we will also encounter bosonic Matsubara frequen- 
cies, which we denote by fl n = lixTn. Using again a four- 
vector notation, k = (k,w„), J^k = T X^„ / dk/(2n) d , 
we can write, 

S* =^^V^(fc) [^„-k 2 /2m + M ] V.(fc) , (2.2a) 



Si n t = — ^2 ^2 (5 fei+fc2,fc 3 +fc4 v ( k 2 - h) 

<7l,&2 {ki} 



(2.2b) 



For the long- wavelength, low-frequency processes we will 
be interested in, only the scattering of particles and holes 
close to the Fermi surface is important. It is custom- 
ary anii-convenient to divide these processes into three 
classesJla (1) Small-angle scattering, (2) large-angle scat- 
tering, and (3) 2fcp-scattering. These classes are also re- 
ferred to as the particle-hole channel for classes (1) and 

(2) , and the particle-particle or Cooper channel for class 

(3) , respectively. The correspondingscattering processes 
are schematically depicted in Fig. |jj For our purposes 




,(2) 




FIG. 1. Typical small-angle (1), large-angle (2), and 
2fcF-scattering processes (3) near the Fermi surface in d = 2. 

it is convenient to make the phase space decomposition 
that is inherent to this classification explicit by writing 
the interaction part of the action, 



.(1) , o( 2 ) i c( 3 ) 
J int ' °int °int 



(2.3a) 



where 



-T 



c(2) _ 



w ( q ) (k)^ (p + q) 

&i,o'2 k,p q 

xip<T2(p)i><ri(k + q) , 
-T 

(71,(72 k,p q 

xip a2 (k + q)ip ai (p) 



(2.3b) 

J ' / _ _ 

T EEE^P - ^ ^ {k)^a 2 (p + q) 



s, 



(3) 



-T 



(2.3c) 

J212' v ( k + V) ^{k)^{-k + q) 

(Ti^ct 2 k,p q 

xiPaAp + qW<n(-p) ■ ( 2 -3d) 



Here the prime on the q-summation indicates that only 
momenta up to some cutoff momentum A are integrated 
over. This restriction is necessary to avoid double coun t- 
i ng, s ince each of the three expressions, Eqs. ( [2.3b ) - 
( |2.3d| ) , represents all of Sj nt if all wavevectors are summed 
over. The long-wavelength physics we are interested in 
will not depend on A. 

The above phase space decomposition is correct in di- 
mensions d > 2. In d = 1, the Fermi surfaces col- 
lapses onto two Fermi points, and the processes we called 
above large-angle scattering and 2fcp--scattering become 
indistinguishable. The three independent scattering pro- 
cesses are usually chosen as the ones shown in Fig. ||, 
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p + q, G 2 



p,G 2 



p + q,0 2 k + q,G 2 



FIG. 2. The three independent scattering processes near 
the Fermi surface with interaction amplitudes gi, 32, and 34 
in d = 1. 



k, a 1 



and the corresponding coupling interaction potentials are 
denoted by g±, g 2 , and 34." Inspection shows that the 
action written in Eqs. ( |2.3| ) counts each of these pro- 

(3) 

cesses twice. If S^t i s dropped, then the 34-process is 
still counted twice. However, it is known that 34 does 
not-contribute to the logarithmic terms we are interested 
mH3 For our purposes it therefore is sufficient to just 
drop the particle-particle channel when we are dealing 
with d = 1. 



p + q,C7 2 -k + q, G 2 

FIG. 3. The three interaction vertices with coupling con- 
stants Ti, T2, and T3. 



B. Simplifications of the Model 

Th e ef fectiv e interaction potentials that appear in Eqs. 
( |2.3b| ) - ( 2.3d| ) are all given by the basic potential v, taken 



at different momenta. S int contains the direct scatter- 
ing contribution, or v(q), with q the restricted momen- 
tum. If v is chosen to be a bare Coulomb interaction, 
then this leads to singularities in perturbation theory in 
v that indicate the need for infinite resummations to in- 
corporate screening. For simplicity, we assume that this 
procedure has already been carried out, and take v to 
be a statically screened Coulomb interaction. For effects 
that arise from small values of |q| it is then sufficient to 
r epla ce v(q) by the number Ti = w(q — > 0)£j In Eqs. 
( 2.3c ), ( 2.3d ) the moduli of k and p are equal to hp for 
the dominant scattering processes, and one usually ex- 
pands these coupling constants in Legendre polynomials 
on the Fermi surface. While all of the terms in this ex- 
pansion contribute to the processes we want to study, we 
note that the coefficients in the angular momentum ex- 
pansion are independent coupling constants. In order to 
establish the existence of a nonanalytic term in Xs(Q)i it 
therefore is sufficient to establish its existence in a par- 
ticular angular momentum channel. For simplicity we 
choose the zero angular momentum channel, I = 0. We 
then have three coupling constants in our theory, namely 
Ti, and T 2 and r 3 , which are v(k — p) and v(k + p), 
respectively, averaged over the Fermi surface. Instead of 
Ti and T 2 one often uses the particle-hole spin singlet 
and spin triplet interaction amplitudes T s and T t that 



are linear combinations of Ti and T 2 . They are related 
to the Fermi liquid parameters Fq and Fq by, 



r s = r a -r 2 /2 = 



1 



2N F 1 + F£ 



v t - r 2 /2 = 



-1 



^0 



2N F 1 + Fq 



(2.4a) 



(2.4b) 



where Np is the density of states at the Fermi level. Our 
simplified model is tantamount to taking only Fq and 
Fq into account instead of the complete sets of Landau 
parameters. As explained above, this is sufficient for our 
purposes. We also define the Cooper channel amplitude, 



r 3 /2 



(2.4c) 



and again we keep only the I = channel. The particle- 
particle channel is neglected in Land au t heory. 



Our model is now defined as Eqs. (2^) and (2^3), with 
^(1)7 V {P " k), and v(k + p) replaced by Ti, T 2 , and 
T3, respectively. We thus have three different interaction 
vertices that are shown in Fig. ||. In the following section 
we will calculate \s in perturbation theory with respect 
to the interaction amplitudes Ti, T 2l and T3. 



III. PERTURBATION THEORY 

A. Contributions to second order in the interaction 

We now proceed to calculate the spin susceptibility, 
Xs, in perturbation theory with respect to the electron- 



(1) 









(6) 



FIG. 4. Second order diagrams that contribute to the non- 
analytic behavior of \ s . The solid vertical line denotes the 
external spin vertex a. 

electron interact|ipju-.This can be done by means of stan- 
dard methodsOElo We will be interested only in con- 
tributions that lead to a nonanalytic wavenumber de- 
pendence. It is easy to see that no nonanalytic behavior 
can occur at first order in the interaction. At second or- 
der, there is also a large number of diagrams for which, 
this is true, and others vanish due to charge neutrality.EHl 
There remain seven topologically different second order 
diagrams, all shown in Fig. 0, that need to be considered. 
We thus write 
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Xs(Q) = 2xo(Q) + y~^X^ + (analytic contributions) , 

i=l 

(3.1) 

where xo denotes the Lindhard function, and the correc- 
tion terms are labeled according to the diagrams in Fig. 
|]. Here and in the remainder of this section we use again 
the four- vector notation of Sec. y, so Q = (Q, Q„), etc. 

These diagrams can be expressed in terms of integrals 
over electronic Green's functions, or b are electron prop- 
agators, that can be read off Eq. (2.2a), 



1 



iLu n — k 2 /2m + /i 
In terms of the Gk , we find 



(3.2) 



X 



(i) 



.-(2) 



J {3) (q,Q) 



(3.3b) 



-2Ti r 2 J2 v 2 J2' ji3) (?> Q) J {3) (-9: -Q) 



(3.3c) 

X (4) = (r 3 ) 2 £ a l a 2 {l-5 rTin )Y^I ( t\q,Q)I^\q,Q), 
&i,<72 q 

(3.3d) 



/(5) 



X-' = (r 3 ) 2 J2 ^(l-^J^I^foQ)!^) , 

CT1,<7 2 q 

(3.3e) 

x (6) -2 *?£' [( r i) 2 J i (4) ( < 7,<3) + (r 2 ) 2 j[ i] {-q,Q) 
xJ (2 \-q) + 2(r 3 ) 2 a & (! - <W=) 



A 7 ) 



J2'4 4) (q,Q)I {2 Hq) , 



(3.3f) 



+(r 2 ) 2 (j(3)(_ g , Q )) ; 



x 



(4) 



(3-3g) 



Here q is a bosonic frequency-momentum integration 
variable. In Eqs. (p. 3]), the following multiplication fac- 
tors have been taken into account. In diagram (1), cither 
one of the interaction lines can be a Ti; the other one 
is then necessarily a T 2 . This leads to a multiplication 
factor of 2, and another factor of 2 comes from the ex- 
istence of an equivalent symmetric diagram. In diagram 
(2), again either one of the two interaction lines can be a 
Ti, with the other line then being a T 2 , but here the two 
expressions one obtains are not identical. Again, there 
is an overall symmetry factor of 2. The same holds for 
diagram (3), but without the overall symmetry factor. 
Diagrams (4) and (5) can be realized only with T3, and 
they carry no multiplication factors. In diagrams (6) and 
(7), both interaction lines must be the same, and diagram 
(6) carries an extra symmetry factor of 2. The spin struc- 
tures represent the fact that the interaction cannot flip 
the spin, and that the external vertex carries a fa ctor 
of a. The functions in the integrands of Eqs. (3J5) are 
defined as, 



J^{q)=Y,GkG 



k — q 



(3.4a) 
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k 

j[ i) {q,Q) = Y J {G k f Gk- q G k -Q 

k 

(l, Q) = Gk Gk-q Gk-Q Gk-q-Q 

k 

I { - 2 \q) = Y;GkG- k+q 
k 

I ( 1 3) (q,Q) = J2 G -« G k+ q G-k-Q 
k 

I^ili Q) = ^ G k G-k+q G -k+q-Q 
k 



lf\q,Q) = Y, G k G-k+q G k+Q G- 



k+q-Q 



(3.4b) 
(3.4c) 
(3.4d) 
(3.4e) 

(3.4f) 
(3-4g) 
(3.4h) 

(3.4i) 



The info rm ation we are interested in is contained in 
Eqs. (3.1) - ( |3.4[ ) in terms of integrals. The remain- 
ing task is to do these integrals. While it is easy to see 
by power counting that all of the above contributions to 
Xs do indeed scale like for 1 < d < 3, and like 

0(1) and 0(Q 2 ) with logarithmic corrections in d = 1 
and d = 3, respectively, we have found it impossible to 
analytically do the integrals in general, i.e. for a finite ex- 
ternal wavenumber in arbitrary dimensions d. However, 
for a perturbative confirmation of the expected nonana- 
lyticity such a general analysis is not necessary. Rather, 
it is sufficient to explicitly obtain the prefactors of the 
logarithmic singularities in d = 1 and d = 3. If they 
are not zero, then by combining this with power count- 
ing and the expected continuity of \s as a function of d, 
it follows that the prefactor of the nonanalyticity 
does not vanish for generic values of d either. For the 
temperature dependence at Q = the integrals can be 
done in arbitrary d, see Sec. 



In Sees. IIIB 
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below. 

HID we therefore analyze the above 
integrals in d = 1 and d = 3. In doing so, we treat the 
particle-hole and particle-particle channel contributions 
separately, since they have quite different structures. We 
also anticipate that we will be interested only in the static 
spin susceptibility, so Q — (0, Q). In d = 1, we write Q 
for the one-dimensional vector, i.e. a real number that 
can be either positive or negative. 



B. d= 1 

Let us first consider d = 1. We do this mainly to make 
contact with established results in the literature. As ex- 
plained above, the particle-particle channel must not be 
taken into account in d = 1, so we put T3 = 0. Since we 
are interested in a logarithm that results from an infrared 
singularity, it suffices to calculate the integrands in the 



limit of small frequencies and wa venum bers. Be perform- 
ing the integrals in Eqs. ( 3.4a ) - ( 3.4d ) one obtains, with 
Q = (0,Q) and q= (fi n ,q), 



J (2) (q) = 



-N f 



1 + (n n /v F qy 



(3.5a) 



J^(q,Q) = N f 



jf ) {q,Q) = N I 



if2 n q/Q 



ifi n (Q-q)/Q 



q (v F q) 



Q (na 
q 2 /Q 2 



(Q 



MQ-q)) 2 . 

(3.5b) 

q/Q 

-q)7Q 2 



4 4) (q : Q)=N F 



2q 2 /Q 2 



(3.5c) 

(g-Q)7Q 2 

W n + (v F (q-QW 
(3.5d) 



Inserting this into Eqs. ( p.3[ ), doing the final integrals, 
and collecting the results one obtains, apart from analytic 
terms, 



ni + (v F qy 

(q+Q)7Q 2 



" 2 + (Mq + Q)) 2 



X.(Q) = 2N F - 4N F (T t N F ) 2 \n(2k F /\Q\) 



(3.6) 



This msult agrees with the well known one to this order 
in r t .Efl One would expect that the In |Q| gets replaced 
by a In O or InT if one works at Q = and finite Q or 
T, respectively. Explicit calculations confirm this. Of 
course the physical content of this perturbative result is. 
limited, since the ground state is not a Fermi liquid.ESI 
For later reference we also mention that, to logarithmic 
accuracy, it is not necessary to keep Q nonzero in the 
above calculation. If one works at Q = and determines 
the prefactor of the resulting logarithmic divergence, then 
one obtains the same result as above. 



C. Particle- hole channel in d = 3 

In d = 3, both the particle- hole and the particle- 
particle channel contribute to the terms we are inter- 
ested in. Since the structures of the integrals in the two 
channels are quite different, we first consider the particle- 
particle channel. In d = 3, the logarithm appears only 
at 0(Q 2 ). Keeping Q explicitly in the integrals to that 
order would be hard. However, as was pointed out in 
the preceding subsection, to logarithmic accuracy this is 
not necessary. Rather, we can just expand in Q. The 
prefactor of the Q 2 term will then be logarithmically di- 
vergent, and the prefactor of the divergence will be the 
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same as that of the Q 2 In | Q | term whose pres enc e the 
divergence signalizes. By expanding Eqs. ( |3.4b| ) - ( 3.4d ) 
to 0(Q 2 ), and dropping the uninteresting contribution 
to the homogeneous Xs, we can express all logarithmic 
contributions to Xs in terms of two integrals, 



(Gk+q) Gk ^ ' Gp Gp-g 

/ p 

• 




D. Particle- particle channel in d — 3 

We now turn our attention to t he p article-particle 
channel. As can be seen from Sec. II A , diagra ms (4 ) 
- (7) in Fig. | contribute. From Eqs. Q3.3dp and Q3.3^ ) 
it follows that the particle-particle channel contributions 
of diagrams (4) and (7) cancel each other, so we are left 



with y^ 5 ) a nd x^ & \ ■ Expanding the functions l[ 4 ^ and 7 L 
Eqs. (3.4h), (|3.4j), to order Q 2 and doing the integrals 



(4) 



one finds that the leading logarithmic contributions to 
both x^ an d X*- 6 -* can be expressed in terms of a single 
integral, 



^lE'E m (G k f G k ^{Gpf G p , 



4^ ^ \ m 

q k 
= ~Jl 



(3.7b) 



where we have kept only the most divergent term. We 
find, 



/(2) 



x (1) = -sr^ q 2 Ji 



-4r 1 r 2 q 2 (Ji + j 2 ) 



c (3) = -sr^ q 2 j 2 



(3.8b) 



(3.8c) 



A 7 ) 



8(r 2 + r 2 ) q 2 Ji 



q 2 (-r 2 Ji - r 2 j 2 ) 



(3.8d) 



(3.8e) 



Here we have used the fact that the structure (J^ 3 -*) 2 that 
appears in x ; X^i an d x' 7 \ if expanded to order Q 2 , 
yields two terms, one of which gets canceled by parts of 
the other. The remaining contribution can be expressed 
in terms of J 2 . 

We see that in the skeleton diagrams, x^ 1 ' - x' 3 \ self 
energy contributions and vertex corrections cancel each 
other. However, in the insertion diagrams, x an d X^\ 
the same cancellation is effective only in the spin singlet 
channel, while in the spin triplet channel the two dia- 
grams add up. Interpreting the logarithmic divergence 
in Ji as a ln|Q| as explained above, we obtain for the 
particle- hole channel contribution to Xs, 



f- h = 2N F + 2N F (T t N F ) 2 - 



_Q 

2k F 



ln(2fc F /|Q|) 



(3.9) 



/ = E'E 



k Q 



q k 



(Gk) G-k+q ^ ' G p G-p+q 
P 

(3.10) 



Inspection of the integrand shows that the leading di- 
vergency in / is a logarithm squared, in contrast to the 
particle-hole channel, where the leading term is a simple 
(3.8a) logarithm. The reason is that G p G- p + q contains a 
term ~ In |q| for q — > 0, which is just the usual BCS-type 
logarithm that is characteristic of the particle-particle 
channel. It also depends on an ultraviolet cutoff, since 
Gp G-p+ q does not exist in d = 3 if the integration is 
extended to infinity. In conjunction with the other factor 
in the integrand of /, which is an algebraic function, this 
gives the leading behavior 



dq ln|q| 



i(uj/v F ) 



o [ q 2 + (uj/v f ) 2 



(3.11) 



While this diverges like (InO) 2 by power counting, the 
prefactor of the divergency turn s out to be zero since the 
frequency integral in Eq. (3.11) vanishes. This leads to 
the following conclusion for the particle-particle channel 
contribution to Xs, 



Xs 



2N F + 2N F (T C N F ) 2 x (ln(2fc F /|QQ) 



+OQn(2k F /\Q\)) 



(3.12) 



Our method of expanding in powers of Q, and extract- 
ing the prefactor of the ensuing singularity, works only for 
the leading nonanalytic contribution. With th is m ethod, 
therefore, the result that is expressed in Eq. ( 3.12 ) is all 
we can achieve. In order to determine the prefactor of 
the next-leading In | Q | term, one would have to keep a 
nonzero external wavenumber explicitly. As pointed out 
before in the context of the particle-particle channel, this 
would be very difficult. However, for our purposes this is 
not really necessary. We know that the interaction ampli- 
tudes in the particle-hole and particle-particle channels, 
respectively, are independent. Therefore, the particle- 
particle channel contribution cannot in general cancel the 
nonzero contribution from the particle-hole channel that 
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we found in Sec. IIIC. What we have established is that 



the particle-particle channel is not more singular than 
the particle-hole channel, and for showing that the lead- 
ing nonanalyticity in Xs is In |Q| with a nonzero prefactor 
this is sufficient. 

It should be pointed out that low-order perturbation 
theory probably overestimates the importance of the 
particle-particle channel. Usually, singularities in the 
particle-particle channel are logarithmically weaker than 
those in the particle-hole channel, since a BCS-type lad- 
der resummation changes a In a; singularity into a minx, 
and a x v singularity into a x y / In x. We expect this mech- 
anism to work in the present problem, so the particle- 
particle channel singularities are probably in fact asymp- 
totically negligible compared to the particle-hole channel 
ones. We also note that so far we have not really estab- 
lished that higher order terms in the perturbation expan- 
sion cannot lead to stronger singularities than the ones 
we found at second order in the interaction amplitudes. 
This point will be further discussed in Section IV below. 



E. Temperature dependence of Xs(Q — 0) 

In the last two subsections we have established that 
Xs in d = 3 at T — does indeed have a nonanalytic 
contribution proportional to Q 2 ln|Q|. As we pointed 
out in the introduction, in a Fermi liquid the wavenum- 
ber scales like frequency or temperature, and one would 
therefore naively expect a T 2 In T contribution to the ho- 
mogeneous Xs at T > 0. This raises the question of 
whether our psults are compatible with those of Carneiro 
and Pethick,E3 who did not find such a contribution. In 
order to clarify this, let us calculate Xs(Q = 0, T) explic- 
itly within our formalism. For the reasons explained in 
Sec. HID we restrict ourselves to the particle-hole chan- 



nel, as did Ref. |13j. 

To this end, we put Q = in Eqs. ( jOel) - ( |3.4d| ), and 
nsider the temneratnre denendence nf vW - v^- 1 . v( 6 ) 



consider the temperature dependence of x y ^' - X y "'i X 
and x^ ■ The relevant integrals are of the structure, 



q 2 T V /(q, in n ) g(q, ifi„ 



iS2„ 



(3.13) 



which are most conveniently done by using the spec- 
tral representation for the causal functions /(q, iO n ) and 
g(q, ifl„)E3 Simple considerations show that there is no 
T 2 In T term if both / and g are algebraic functions; only 
if at least one of them possesses a branch cut can such 
a nonanalyticity arise. This immediately rules out x^\ 
and the first and second contribution to x and X^\ re ~ 
spectively, as sources for a T 2 In T. The rea son is that an 
explicit calculation of J^\q,Q — 0), Eq. ( 3.4b| ), in the 
limit of small q shows that the only singularities in this 
function are poles. The same is true for j[ A) (q,Q = 0) 

and 4 4) (q,Q = 0), but J (2) (g), which is minus the Lind- 
hard function, has a branch cut, and so all of the remain- 
ing terms potentially go like T 2 lnT. 



Since again we are aiming only at logarithmic accu- 
racy, we can replace (q, Q = 0) and 4 (q, Q = 0) 
by low-frequency, long-wavelength expressions for which 
4 4) (<1'Q = 0) = -2J 1 (4) (g,g = 0). The contributions 
from x^ and x therefore cancel (remember that dia- 
grams (1) and (2) in Fig. ^ carry multiplication factors 
4 and 2, respectively). The contributions from x^ and 
X^ can both be expressed in terms of an integral 

J= fdWTj24 4 \l,Q = 0)J {2 \q) 



(3.14) 



In doing this integral one may encounter individual terms 
that go like T 2 In T, but all of those terms cancel, and 
the leading T-dependence of J is T 2 . There hence is no 
T 2 InT contribution to Xs in d = 3. 

This result agrees with the conclusion of Ref. [l3| which 
reached it on the basis of Fermi liquid theory. We dis- 
agree, however, with the assertion of that reference that 
within the framework of microscopic perturbation theory 
the absence of the T 2 In T is due to cancellations between 
vertex corrections and self energies, and is hence a con- 
sequence of gauge invariance. What we find instead is 
that, for all diagrams in Fig. f|, the T 2 InT terms vanish 
individually. This is consistent with the result of Ref. |24|. 
These authors calculated Xs in paramagnon approxima- 
tion, which in our language corresponds to taking only 
X^ 6 ' and x^ i n to account, plus infinite resummations 
that contribute to higher orders in the interaction am- 
plitudes. They reported the absence of T 2 lnT terms in 
their calculation, rather than their cancellation between 
the two diagrams. 

This absence of the expected nonanalytic In- 
dependence in d = 3 is somewhat accidental. This can 
be seen from the 1-d case, where, as was pointed out in 
Sec. [II B, there is a InT contribution to the homoge- 



neous spin susceptibility. The technical reason is that in 
d = 1, integrands whose only singularities are poles do 
contribute to the T 2 lnT terms. Consequently, in d = 1 
T and Q are interchangeable in the logarithmic terms, 
while in d = 3 they are not. Furthermore, the same type 
of integrals that lead to a InT in d = 1, also contribute to 
a T' 1 ^ 1 nonanalyticity in 1 < d < 3. In these dimensions 
we therefore expect to find, 



X 



p -' l (Q = 0) = 2N F + 2N F (T t N F ) 2 c d (T/4e F ) 



d-l 



(3.15) 



with Cd a d-dependent, positive number. 

We also mention that the absence of a T 2 In T in the 
self energy diagrams in d = 3 does not contradict the 
presence of such a term in the specific heat coefficient. 
The relation betweezuthe specific heat and the Green's 
function is intricate^ 2 ! and the resulting integrals have a 
different structure from the ones that determine Xs- 
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IV. DISCUSSION 

A. Our results in a mode-mode coupling theory 
context 

In this section we have a more detailed look at the 
mode-mode coupling arguments that were presented in 
the introduction. We also stress some analogies between 
classical and quantum fluids, and discuss some important 
differences between clean and disordered systems. 

Let us consider four distinct systems: (1) a classi- 
cal Lorentz model (i.e. a classical patticle moving in a 
spatially random array of scatterersES) , (2) a classical 
fluid, (3) a Fermi liquid with static impurities, and (4) 
a clean Fermi liquid. These systems represent classical 
and quantum fluids with and without quenched disorder, 
respectively. As was pointed out in the introduction, dy- 
namical correlations are ultimately responsible for all of 
the effects discussed in this paper. However, in classical 
systems they do not manifest themselves in static equi- 
librium properties, while in quantum systems they do. 
In order to discuss the analogies between classical and 
quantum systems, let us therefore digress and consider 
an equilibrium time correlation function. A convenient 
choice is the current-current correlation function, whose 
Fourier transform determines the frequency dependent 
diffusivity D(il). In both of the classical systems, (1) and 
(2), this correlation function exhibits a long-time tail, so 
D(O) is nonanalytic at = 0. For O — > one finds for 
the classical Lorentz model, 



k+q, CO + Q. 



D(n) /d(o) = i + a in + b {in) d/2 



(4.1a) 



while for the classical real fluid one finds, 

D(fi)/D(0) = 1 - b' {in) id - 2 ^ 2 . (4.1b) 



The coefficients b and b' in Eqs. (LI) are positive. The 
long-time tail in the real fluid is stronger than the one in 
the Lorentz gas because the former has more soft modes. 
More importantly, the static scatterers in the Lorentz gas 
lead to a sign of the effect that is different from the one 
in the real fluid. All of these features can be understood 
in terms of the number and the nature of the soft modes 
in these systems.cl In disordered Fermi liquidsu one has, 



D(Q)/D(0) = 1 + 6" (in) {d - 2)/2 



(4.1c) 



with b" > 0. Here the sign is the same as in the classical 
Lorentz model, which is due to the quenched disorder in 
either system. The strength of the long-time tail, how- 
ever, is equal to that in the classical real fluid. As men- 
tioned in Sec. |[ the coupling of statics and dynamics in 
quantum statistical mechanics leads to a related nonan- 
alyticity in the static spin susceptibility of a disordered 
Fermi liquid, namely 



x s (Q)/x s (0) = i-c|Q| d - 



q,£2 




-k, —CO 

FIG. 5. Mode-mode coupling process describing the decay 
of a current mode (dashed line) into two sound modes (solid 
lines) . 

with c > 0. 

On the basis of these results, it is possible to predict 
both the strength of the singularity, and the sign of the 
prefactor, in the Q-dependence of Xs in a clean Fermi 
liquid, which is what we are mainly concerned with in 
this paper. In order to do so, let us recall the origin of 
the nonanalyticity in the classical fluid, Eq. (4.1b). The 
density excitation spectrum, i.e. the dynamical struc- 
ture factor as measured in a light scattering experiment, 
in a classcical fluid consists of three main features: The 
Brillouin peaks that describe emission and absorption 
of sound waves, and the Rayleigh peak that describes 
heat diffusion. For our purposes, we focus on the for- 
mer. In the density-density Kubo correlation function, 
C(k, lo) (whose spectrum is in a classical system simply 
proportional to the structure factor S*(k, uj)), they man- 
ifest themselves as simple poles J3 



C(k,oj) 



1 



1 



uj — vk + ijk 2 /2 uj + vk — ijk 2 /2 

= C+(k,w) +C_(k,w) , (4.3) 



where v is the speed of sound, and 7 is the sound atten- 
uation constant. Now let us consider the simplest possi- 
blc m ode-mode coupling process that contributes to Eq. 
( 4.1b] ), namely one where a current mode decays into two 
sound modes that later recombine, see Fig. Consider 
a process where one of the internal sound propagators is 
a C+, and the other a C_. At zero external wavenumber, 
this leads to a convolution integral, 

/ duj J dkC+(k,w) CL(-k,-w + fi) 

dk „ 1 . „ ~ r>( rf - 2 )/ 2 . (4.4a) 



1 

n + i'yk 2 

Note that by this mechanism the long-time tail in a sys- 
tem whose low-lying modes have a linear dispersion be- 
comes as strong as the one in a system with diffusive 
modes. In contrast, if both of the sound propagators are 
C4. or C_, one obtains a weaker singularity, 



(4.2) J dflj dk C+(k,uj) C+(-k, 



n) 







dk 



1 



fl - 2vk + iO 



(4.4b) 



Now let us consider the corresponding quantum system, 
i.e. the clean Fermi liquid. Again, the low-lying modes 
(i.e., particle- hole excitations), have a linear dispersion. 
However, at zero temperature the structure factor and 
the Kubo function are no longer proportional to one an- 
other. Rather, the fluctuation dissipation theorem shows 
that they are related by a Bose distribution function that 
eliminates the pole at to = ck from the structure factor. 
This is simply a consequence of the fact that at zero 
temperature there are no excitations that could get de- 
stroyed in a scatter ing p rocess. Consequently, the process 
described by Eq. ( 1.4a ) is not available in this s ystem , 
and one is left with the weaker singularity of Eq. (4.41). 
Since the diffusion coefficient is infinite at T = in a 
clean system, we look instead at the spin susceptibility 
as a function of Q. Q scales like so we expect a singu- 
larity of the form |Q| d_1 , as opp osed to the |Q| d ~ 2 in a 



disordered Fermi liquid, Eq. fl4.2| ). The sign of the prefac- 
tor is determined by whether or not the system contains 
quenched disorder. It should therefore be opposite to the 
sign in the dirty case. We thus expect for the wavenum- 
ber dependence of the spin susceptibility in a clean Fermi 
liquid, 



X s (Q)/Xs(0) = l + c'|Q 



d-l 



(4.5) 



with c' > 0. This is precisely what we found in Sec. Ill 



by means of perturbation theory. Notice that the mode- 
mode coupling arguments suggest that the sign of the 
prefactor c' will be positive, regardless of the interaction 
strength, as is the sign of the long-time tail i n a cl assical 
fluid. We will come back to this point in Sec. IV C below. 



B. Our results in a renormalization group context 

Another useful way to look at our results is from a 
renormalization group point of view. The Fermi liquid 
ground state of interacting fermion systems in d > 1 
has recently been identified with a stable fixed point in 
renormalization group treatments of both a basic fppiion 
theory,E3 and a bosonized version of that theorycJ The 
instability of the Fermi liquid in d — 1 is reflected by 
an infinite number of marginal operators whose scale di- 
mensions are proportional to d — 1, i.e. they all become 
relevant in d < 1, and are irrelevant in d > 1. In the 
present context, the Fermi liquid nature of the ground 
state in d > 1 is reflected by the fact that the homoge- 
neous spin susceptibility is finite in perturbation theory. 
The nonanalytic corrections at finite wavenumber that we 
are interested in correspond to the leading correction to 
scaling in the vicinity of the Fermi liquid fixed point, i.e. 
to an irrelevant operator with respect to that fixed point. 
Among the irrelevant operators, there thus must be one 
whose scale dimension determines the leading wavenum- 
ber dependence of the spin susceptibility. 



An identification of this operator within the framework 
of a renormalization group analysis would not only pro- 
vide another derivation of our result, but would also es- 
tablish that the behavior we have found in perturbation 
theory constitutes the leading Q-dependence to all or- 
ders in the interaction amplitudes. This program has 
not been carried—out yet, although preliminary results 
are encour aging. EJ This will provide a connection be- 
tween the mode-mode coupling arguments presented in 
the previous subsection and renormalization group argu- 
ments that will be analogous to a corresponding connec- 
tion in classical fluids that has been known to exist for 
some time.Q 

In this context it should also be mentioned that there is 
no universal agreement that the ground state of a weakly 
interacting Fermi system in d > 1 is a Fermi liquid. It 
has been proposed that there exists a relevant operator 
that makes the Fermi liquid fixed point .mistable, and 
leads to a non-Fermi liquid ground state.EZI In order to 
destroy the Fermi liquid in d dimensions, this would re- 
quire a long-range effective interaction that falls off more 
slowly than l/r d at large distances. While we do find 
an effective long-range interaction between the spin de- 
grees of freedom, it falls off like l/r 2d_1 , and hence leaves 
the Fermi liquid fixed point intact. The same conclusion 
was reached in Ref. |ll] from studying the specific heat in 
d = 2. 



C. Summary, and physical consequences of our result 

We finally turn to a summary of our results, and to a 
discussion of their physical consequences. By means of 
explicit pcrturbative calculations to second order in the 
interaction, we have found that the wavenumber depen- 
dent spin susceptibility in d = 3 has the form, 



Xs(Q)/Xs(Q = 0) = 1 + c 3 (Q/2fc F ) 2 ln(2fc F /|Q| 



-0(Q 2 ) 



(4.6a) 



We have calculated the particle-hole channel contribu- 
tion to the constant C3, and have found it to be posi- 
tive. More generally, it follows from our analysis that 
in (i-dimensional systems, the spin susceptibility has a 
nonanalyticity of the form 

Xs(Q)/Xs(Q = 0) = l + c d (|Q|/2fc F ) d_1 + 0(Q 2 ) , 

(4.6b) 

where the particle-hole channel contribution to Cd is again 
positive. 

A very remarkable feature of Eqs. (JO]) is the sign of 
the leading Q-dependence: For d < 3, Xs increases with 
increasing |Q| like |Q| d_1 . For any physical system for 
which this were the true asymptotic behavior at small 
Q, this would have remarkable consequences for the zero- 
temperature phase transition from the paramagnetic to 
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the ferromagnetic state as a function of the exchange cou- 
pling. One possibility is that the ground state of the sys- 
tem will not be ferromagnetic, irrespective of the strength 
of the spin triplet interaction, since the functional form of 
Xs leads to the instability of any homogeneously magne- 
tized ground stateBa Instead, with increasing interaction 
strength, the system would undergo a transition from a 
paramagnetic Fermi liquid to some other type of mag- 
netically ordered state, most likely a spin density wave. 
While there seems to be no observational evidence for 
this, let us point out that in d = 3 the effect is only loga- 
rithmic, and would hence manifest itself only as a phase 
transition at exponentially small temperatures, and ex- 
ponentially large length scales, that might well be un- 
observable. For d < 2, on the other hand, there is no 
long-range Heisenberg ferromagnetic order at finite tem- 
peratures, and the suggestion seems less exotic. Further- 
more, any finite concentration of quenched impurities will 
reverse the sign of the leading nonanalyticity, and thus 
make a ferromagnetic ground state possible again. 

Another possibility is that the zero-temperature 
paramagnet-to-ferromagnet transition is of first order. It 
has been shown in Ref . O that the nonanalyticity in \s (Q) 
leads to a similar nonanalyticity in the magnetic equation 
of state, which takes the form 



might lead to a different sign. At this point, neither one 
of these possibilities can be ruled out mathematically. 
However, from a physical point of view neither one is 



tm — Vdin + u m 3 = h 



(4.7) 



with m the magnetization, h the external magnetic field, 
and «>0a positive coefficient. If the soft mode mech- 
anism discussed above is the only mechanism that leads 
to nonanalyticities, then the sign of the remaining coeffi- 
cient v in Eq. ( 4.7 ) should be the same as that of Cd in Eq. 
( 4.6b| ), i.e. Vd > 0. This would imply a first order tran- 
sition for 1 < d < 3. In this case the length scale that in 
the previous paragraph would have been attributed to a 
spin density wave would instead be related to the critical 
radius for nucleation at the first order phase transition. 
Further work will be necessary to decide between these 
possibilities. 

The conclusion that there is no continuous zero tem- 
perature paramagnet-to-ferromagnet transition is in- 
escapable for any system with a particle-hole channel in- 
teraction that is sufficiently weak for our perturbative 
treatment to be directly applicable. An important ques- 
tion is now whether or not it holds more generally, for sys- 
tems whose interactions are in general not weak. There 
are four obvious mechanism by which the sign of the lead- 
ing Q-dependence of Xs could be switched from positive 
to negative: (1) Higher order contributions could lead to 
a sign of Cd for realistic interaction strengths that is differ- 
ent from the one for weak interactions, or (2) they might 
lead to a stronger singularity with a negative prefactor 
that constitutes the true long-wavelength asymptotic be- 
havior, or (3) the particle-particle channel contribution 
might have a negative sign that overcompensates the pos- 
itive contribution from the particle-hole channel, or (4) 
the higher angular momentum channels that we neglected 



very likely to occur. As we have explained in Sec. IV A 



both the functional form and the sign of the nonanalyt- 
icity found in perturbation theory are in agreement with 
what one would expect on the basis of a suggestive anal- 
ogy with classical fluids. Also, the renormalization group 
arguments sketched in Sec. IV B make it appear likely 



that Eqs. (4.6) constitute the actual asymptotic small- 
Q behavior of Xsi although an actual renormalization 
group proof of this is still missing. This makes the first 
two possibilities appear unlikely. The third possibility is 
unappealing for two reasons. First, the effective inter- 
action in the particle-particle channel is typically much 
weaker than the one in the particle-hole channel. The 
reason is the characteristic ladder resummation that oc- 
curs in the particle-particle channel if one goes to higher 
orders in perturbation theory. This leads to an effective 
interaction of the 'Coulomb pseudopotential' type that 
is much weaker (typically by a factor of 5 to 10) than, 
what low order perturbation theory seems to suggest .E3 
Second, that same resummation weakens any singular- 
ity (cf. the discussion at the end of Sec. HID), which 



probably makes the particle-particle channel singularity 
subleading. Finally, the higher angular momentum Fermi 
liquid parameters are usually substantially smaller than 
the ones at I = 0, which makes possibility (4) unlikely, 
except possibly in particular systems. 

If the sign of the nonanalyticity is, for some reason, 
negative at the coupling strength necessary for a ferro- 
magnetic transition to occur, at least in some systems, 
then in these systems the quantum phase transition from 
a paramagnet to a ferromagnet at zero temperature as a 
function of the interaction strength will be a conventional 
continuous quantum phase transition with an interesting 
critical behavior. This is because the nonanalyticity in 
Xs leads to an effective long-range interaction between 
spin fluctuations, which in turn leads to critical behavior 
that is not mean-field like, yet exactly solvable. This has 
been discussed recently in some detail!! 
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